Lennard— Jones and Lattice Models of Driven Fluids 
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We introduce a nonequilibrium off-lattice model for anisotropic phenomena in fluids. This is a 
Lennard- Jones generalization of the driven lattice-gas model in which the particles' spatial coordi- 
nates vary continuously. A comparison between the two models allows us to discuss some exceptional, 
hardly realistic features of the original discrete system — which has been considered a prototype for 
nonequilibrium anisotropic phase transitions. We thus help to clarify open issues, and discuss on 
the implications of our observations for future investigation of anisotropic phase transitions. 
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The concept of nonequilibrium phase transition (NPT) 
IH Q< 01 helps our present understanding of many com- 
plex phenomena including, for instance, the jamming in 
traffic flow on highways [4J , the origin of life [5| , and the 
pre-humans transition to mammals [fjl. Many studies of 
NPTs have focused on lattice systems|lllla[l3,[ll|,[ll 
Il3l ]. This is because lattice realizations are simpler than 
in continuum space, e.g., they sometimes allow for exact 
results and are easier to be implemented in a computer. 
Furthermore, a bunch of emerging techniques may now 
be applied to lattice systems, including nonequilibrium 
statistical field theory. A general amazing result from 
these studies is that lattice models often capture the es- 
sentials of social organisms, epidemics, glasses, electrical 
circuits, transport, hydrodynamics, colloids, networks, 
and markets, for example. 

The driven lattice gas (DLG) [T^. a crude model 
of "super ionic currents" [jjj , has become the theoreti- 
cians' prototype for anisotropic NPTs. The DLG con- 
sists of a lattice gas with the particles hopping prefer- 
entially along one of the lattice directions, say x. One 
may imagine this is induced by an external drive, Ex, 
e.g., an applied electric field assuming the particles are 
positive ions. Consequently, for periodic boundary con- 
ditions, a particle current and an anisotropic interface 
set up along x at low temperature, T < Te- That is, 
a liquid-like phase which is striped then coexists with 
its gas. More specifically, assuming — for simplicity and 
concreteness — the square lattice half filled of particles, 
Monte Carlo (MC) simulations show that the function 
Te monotonically increases with E from the Onsager 

value T = Ibnsager = 2.269Jfc^ 1 to T^ ~ 1.41bnsager- 

This limit corresponds to a nonequilibrium critical point. 
As a matter of fact, it was numerically shown to belong 
to a universality class other than the Onsager one, e.g., 
MC data indicates [3 ~ 0.33 (instead of the Ising value 
P = 1/8 in two dimensions) for the order parameter crit- 
ical exponent [U El 13 El 

Statistical field theory is a complementary approach 
to the understanding of nonequilibrium ordering in the 
DLG. The derivation of a general mesoscopic description 



is still an open issue, however. Two different approaches 
have been proposed. The driven diffusive system (DDS) 
H H^. |20| , which is a Langevin type of equation aimed 
at capturing all the relevant symmetries, predicts that 
the current will induce a predominant mean-field behav- 
ior and, in particular, /3 = 1/2. The anisotropic diffusive 
system (ADS) 21], which follows after a non-rigorous 
coarse graining of the master equation, rules out the rel- 
evance of the current and leads to the above-indicated 
MC critical exponent for E — > 00. However, the ADS ap- 
proach reduces to the DDS for finite E, a fact which is 
hard to be fitted to MC data, and both contain disquiet- 
ing features [22| . 

In any case, field theoretical studies have constantly 
demanded further numerical efforts, and the DLG is 
nowadays the most thoroughly studied system showing 
an anisotropic NPT. The topic is not exhausted, how- 
ever. On the contrary, there remain unresolved matters 
such as the above mentioned issues concerning critical 
and mesoscopic behaviors, and the fact that Te increases 
with E, which is counterintuitive |23l ]. Another signi- 
ficative question concerns the observation of triangular 
anisotropics at early times after a rapid MC quench from 
the homogenous state to T < Tg. The triangles happen 
to point against the field, which is contrary to the pre- 
diction from the DDS continuum equation [24j. 

This paper reports on a new effort towards better un- 
derstanding basic features of NPTs. With this aim, we 
here present a description of driven systems with con- 
tinuous variation of the particles' spatial coordinates — 
instead of the discrete variations in the DLG. We hope 
this will provide a more realistic model for computer sim- 
ulation of anisotropic fluids. Our strategy to set up the 
model is to follow as closely as possible the DLG. That 
is, we analyze an off-lattice representation of the DLG, 
namely, a microscopically continuum with the same sym- 
metries and hopefully criticality. Investigating these 
questions happens to clarify the puzzling situation indi- 
cated above concerning the outstanding behavior of the 
DLG. The clue is that the DLG is, in a sense, patholog- 
ical. 
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FIG. 1: Typical configurations during the stationary regime 
of the DLJF subject to a horizontal field of intensity E* = 1. 
These graphs, which are for N = 900 particles and density 
p* = 0.30, illustrate, from left to right, (i) coexistence of a 
solid and its vapor (the configuration shown is for tempera- 
ture T* — 0.20), (ii) liquid-vapor coexistence (T* = 0.35), 
and (Hi) a disordered, fluid phase (T* = 0.50). The left-most 
graph shows a detail of the solid strip. The particles, which 
move in a square box of side y N/ p* , are given here an arbi- 
trary size, namely, diameter= l.lcr. 
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FIG. 2: Schematic phase diagrams for the three models, 
as defined in the main text, showing ordered (OP) and dis- 
ordered (DP) phases. The left graph is for the DLG, for 
which To = Tbnsagor- The graph on the right is valid for both 
the NDLG, i.e., the DLG with next-nearest-neighbor (NNN) 



hops, for which To 
To = 0.459 [H . 
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and the DLJF, for which 



Consider a fluid consisting of N particles in a two- 
dimensional L x L box with periodic boundary condi- 
tions. Interactions are according to a truncated and 
shifted Lennard-Jones (LJ) potential [25| : 
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Here, = \fl — fj\ is the relative distance between par- 
ticles % and j, Vlj(t) = 4e [(a/r) 12 — (cr/r) 6 ] , e and fl- 
are our energy and length units, respectively, and r c is 
the cut-off that we shall fix at r c = 2.5a. The preferen- 
tial hopping will be implemented as in the lattice, i.e., by 
adding a drive to the potential energy. Consequently, the 
familiar energy balance is (assuming ks = 1 hereafter): 



—T~ 



(H(?)-H(c)) + Ex-8 



(1) 



Where c = {f±, ...,r]v} stands for configurations, H(c) = 
'52i < jV(rij), and S = rj — f\ is the attempted parti- 
cle displacement. Defining the latter is a critical step 
because, as we shall show, the resulting (nonequilib- 
rium) steady state will depend, even qualitatively on this 
choice. Lacking a lattice, the field is the only source of 



DLG 



DLJF 




FIG. 3: Triangular anisotropics as observed at early times 
for E — 1 in computer simulations of the lattice (left) and the 
off-lattice (right) models defined in the main text. The DLG 
configuration is for t = 6 x 10 4 MCS in a 128x128 lattice 



with N = 7372 particles and T = 0.4T On 



The DLJF 



configuration is for t = 1.5 x 10 
p* = 0.20 and T* = 0.23. 



MCS, N = 10 4 particles, 



anisotropy, and any trial move should only be constrained 
by a maximum displacement in the radial direction. That 
is, we take < \5\ < A, where A = 0.5a in the simula- 
tions reported here. The temperature T, number density 
p = N/L 2 , and field variables will be reduced according 
to T* = T/e, p* = pa 2 , and E* — Ea/e, respectively. 
Our model thus reduces for E — > to the truncated and 
shifted LJ fluid, one of the most studied models in the 
computer simulation of fluids pH |2(j . 

We studied this driven LJ fluid (DLJF) in the com- 
puter by the MC method using a "canonical ensemble" , 
namely, fixed values for TV, p* , T* , and E* . Simulations 
involved up to N = 10 4 particles with parameters rang- 
ing as follows: 0.5 < E* < 1.5, 0.20 < p* < 0.60, 
0.15 < T* < 0.55. The typical configurations one ob- 
serves are illustrated in Fig. ^ As its equilibrium coun- 
terpart, the DLJF exhibits three different phases (at 
least): vapor, liquid, and solid (sort of close-packing 
phase; see the left-most graph in Fig. QJ. At intermedi- 
ate densities and low enough T* , vapor and a condensed 
phase segregate from each other. The condensed droplet 
(see Fig. ^) is not near circular as it generally occurs in 
equilibrium, but strip-like extending along the field di- 
rection. A detailed study of each of these phases will be 
reported elsewhere [2^ |: we here focus on more general 
features. 

A main observation is that the DLJF closely resembles 
the DLG in that both depict a particle current and the 
corresponding anisotropic interface. However, they dif- 
fer in an essential feature, as illustrated by Fig. |2J That 
is, contrary to the DLG, for which Tg increases with 
E, the DLJF shows a transition temperature Tg which 
decreases with increasing E. The latter behavior was ex- 
pectable. In fact, as E is increased, the effect of the 
potential energy in the balance Eq. becomes weaker 
and, consequently, the cohesive forces between particles 
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FIG. 4: Schematic comparison of the accessible sites a par- 
ticle (at the center, marked with a dot) has for nearest- 
neighbors (NN) and next-nearest-neighbors (NNN) hops at 
equilibrium (left) and in the presence of a large horizontal 
field (right). The particle-hole exchange between neighbors 
may be forbidden (x), depend only on the potential energy 
(T), or occur with probability 1 (E). 



tend to become negligible. Therefore, unlike for the DLG, 
there is no phase transition for a large enough field, and 
Te — > for E — > oo in the DLJF. Confirming this, typi- 
cal configurations in this case are fully homogeneous for 
any T under a sufficiently large field E. One may think 
of variations of the DLJF for which Te—,<x> = const > 0, 
which is more realistic, but the present one follows more 
closely the DLG microscopic strategy based on Eq. 

Concerning the early process of kinetic ordering, one 
observes triangular anisotropies in the DLJF that point 
along the field direction. That is, the early-time 
anisotropies in the off-lattice case (right graph in Fig. |3J) 
are similar to the ones predicted by the DDS, and so they 
point along the field, contrary to the ones observed in the 
discrete DLG (left graph in Fig.[3J). 

The above observations altogether suggest a unique ex- 
ceptionality of the DLG behavior. This is to be associ- 
ated with the fact that a driven particle is geometrically 
restrained in the DLG. In order to show this, we stud- 
ied the lattice with an infinite drive extending the hop - 
ping to next-nearest-neighbors (see also Refs.[28l |29|). 
As illustrated in Fig. this introduces further relevant 
directions in the lattice, so that the resulting model, 
to be named here NDLG, is expected to behave closer 
to the DLJF. This is confirmed. For example, one ob- 
serves in the discrete NDLG that, as in the continuum 
DLJF, Te decreases with increasing E — though from 
To = 2.35Tonsagcr in this case. This is illustrated in 
Fig.0 

There is also interesting information in the two-point 
correlation function and its Fourier transform, S(k). In 
fact, the DLG displays (more clearly above criticality) 
slow decay of two-point correlations 30] due to the de- 
tailed balance violation for E ^ 0. For a half-filled lat- 
tice, this function is C(r) = (sgs?±g) — 1/4, where s$ 
is the occupation number at site q and the steady av- 
erage (...) involves also averaging over q. Analysis of 
the components along the field, C(x,0), and transverse 
to it, C(0, y), shows that correlations are qualitatively 
similar for the DLG and the NDLG — although some- 
what weaker along the field for NNN hops. That is, 
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FIG. 5: Parallel (circles) and transverse (triangles) compo- 
nents of the structure factor above criticality for the DLG on 
a 128 x 128 lattice with NN interactions (filled symbols, left 
scale) at T = 1.58 and NNN interactions (empty symbols, 
right scale) at T = 1. Temperatures normalized to Tbnsager. 
Here, n x , y = 12&k x , y /2-k . 



allowing for a particle to surpass its forward neighbor 
does not modify correlations. The power-law behav- 
ior translates into a discontinuity of S{k) [3(j, namely, 
limfe^^o S(k x , 0) 5^ ]imk y —>o S(0, k y ), which is clearly con- 
firmed in Fig.0for both NN and NNN hops. 

The above shows that the nature of correlations is not 
enough to determine the phase diagram. There are al- 
ready indications of this from the study of equilibrium 
systems, and also from other nonequilibrium models. 
That is, one may have a qualitatively different phase di- 
agram but essentially the same two-point correlations 
by modifying the microscopic dynamics [53, . It fol- 
lows, in particular, that the exceptional behavior of the 
DLG cannot be understood just by invoking the func- 
tions C(r) and S(k) or crude arguments concerning sym- 
metries. The fact that particles are constrained to travel 
precisely along the two principal lattice directions in the 
DLG is the cause for its singular behavior. Allowing 
jumping along intermediate directions, as in both the 
NDLG and the DLJF, modifies essentially the phase dia- 
gram but not features — such as power-law correlations — 
that seem intrinsic of the nonequilibrium nature of the 
phenomenon. 

Rutenberg and Yeung [2^ also performed quenching 
experiments for variations of the DLG. They showed, in 
particular, that minor modifications in the DLG dynam- 
ics may lead to an inversion of the triangular anisotropies 
during the formation of clusters which finally condense 
into strips. Our observations above prove that such 
nonuniversal behavior goes beyond kinetics, namely, it 
also applies to the stationary state. 

The unique exceptionality of the DLG has some im- 
portant consequences. One is that this model involves 
features that are not frequent in nature. There are situ- 
ations in which a drive induces stripes but not necessar- 
ily DLG behavior. The fact that a particle impedes the 
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freedom of the one behind to move along x for a large 
enough field may only occur very seldom in cooperative 
transport. The NDLG is more realistic in this sense. In 
any case, the great effort devoted to the DLG during 
two decades has revealed important properties of both 
nonequilibrium steady states and anisotropic phase tran- 
sitions, and there are some unresolved issues yet. For ex- 
ample, a general mesoscopic description which captures 
the exceptionality of the DLG remains elusive. Accord- 
ing to our observations above, such a description needs 
to include the microscopic details of transverse dynam- 
ics which, in particular, should allow one to distinguish 
between the DLG and the NDLG. On the other hand, 
it ensues that, due to its uniqueness, the DLG does not 
have a simple off-lattice analog. This is because the or- 
dering agent in the DLG is more the lattice geometry 
than the field itself. The fact that one needs to be very 
careful when modeling nonequilibrium phenomena — one 



may induce both a wrong critical behavior and an spuri- 
ous phase diagram — ensues again in this example. This 
seems not to be so dramatic in equilibrium where, for 
example, the lattice gas is a useful oversimplification of 
a LJ fluid. 

Finally, we remark that a novel fluid model in which 
the particles move in a continuum space has been intro- 
duced in this paper. The particle 'infinite freedom' is 
realistic, as it is also the LJ potential. It may contain 
some of the essential physics in a class of nonequilibrium 
anisotropic phenomena and phase transitions. On the 
other hand, the model is simple enough to be useful in 
computer simulations, and it is endowed of even-simpler 
and functional lattice analogs such as the NDLG. 

We acknowledge very useful discussions with E. Al- 
bano, M. A. Munoz, and F. de los Santos, and finan- 
cial support from MEyC and FEDER (project FIS2005- 
00791). 
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